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Abstract 



A semi-empirical calculation of the 6s - 7s Stark amplitude a in Cs has 
been performed using the most accurate measurements and calculations of the 
electromagnetic amplitudes available. This is then used to extract the pa- 
rameters of the electroweak theory from experimental data. The results are: 
Q = 269.0(1.3)ao, weak charge of Cs Qw = — 72.41(25)exp(80)t/ieor) deviation 
from the Standard model S = — 1.0(.3)exp(l-0)fheor and limit on the mass of 
the extra Z-boson in 50(10) model Mz^ > 550GeV. 
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Experiments suggested in |I| for measuring parity nonconservation (PNC) in heavy atoms 
have provided an important confirmation f^-Q of the standard model of elementary particles. 
Combining the recent very accurate measurements of parity nonconservation in the Cs atom 
Q with theoretical calculations gives one a possibility to study new physics beyond the 
standard model. The measured nuclear spin-independent part of the PNC effect in Cs is 
of the form (we use the analysis from ) 

- = 1.5939(56)^ (1) 

p cm 

where Ep^c is the PNC El-amplitude of the 6s - 7s transition and j3 is the vector polarizability 
of the transition. The theoretical values of Epj^c 

follows: 

_ n. Qh^ J 0.908(10) Ref [5] 

EpNC - -Ae\aolO (-^) j o.905(9) Ref i 

Here Qw is the weak charge of the cesium nucleus and is the number of neutrons. 

The method for ab initio calculations of Ep^c that we used in ||^ was based on an all- 
orders summation of the dominating diagrams of the many-body perturbation theory in the 
residual Coulomb interaction using a relativistic Hartree-Fock basis set and Green's functions. 
This technique has been described in 

We took into account direct and exchange polarization of the atomic core by the external 
electric field and the weak nuclear potential using the Time-Dependent Hartree-Fock method 
(summation of the "RPA with exchange" chain of diagrams), calculated second-order corre- 
lation corrections and 3 series of dominating higher-order diagrams: 

1. Screening of the electron-electron interaction. This is a collective phenomenon and so 
the corresponding chain of diagrams is enhanced by a factor approximately equal to the 
number of electrons in the external closed subshell (the 5p electrons in Cs). We stress 
that our approach takes into account screening diagrams with double, triple and higher 
core electron excitations §], in contrast to popular pair equations (coupled cluster) 
method, where only double excitations were considered. 

2. Hole-particle interaction. This effect is enhanced by the large zero-multipolarity diagonal 
matrix elements of the Coulomb interaction. 

3. Iterations of the self-energy operator ( "correlation potential"). This chain of diagrams 
describes the nonlinear effects of the correlation potential and is enhanced by the small 
denominator, which is the energy for the excitation of an external electron (in comparison 
with the excitation energy of a core electron). 

The error in the theoretical value was tested in many different ways: by estimating the 
contribution of the unaccounted higher-order diagrams, by comparing the calculated and mea- 
sured values of the energy levels, the fine and hyperfine structure intervals, the probabilities of 
electromagnetic transitions, etc. (see Ref. 0). The result for the PNC amplitude practically 
did not change when we introduced factors into the correlation potential to fit the energy lev- 
els (in imitation of the unaccounted higher-order diagrams). Important tests of our method 
included predictions of the spectrum [|lOl and electromagnetic transition amplitudes for the 



Fr atom |TT|, which is an analogue of Cs. Recently the positions of many energy levels [|T2 



and some transition rates ||T^ of Fr were measured and found to be in excellent agreement 
with our predictions. 

Our calculations of PNC for atoms with electron structures more complex than those of 
the alkaline atoms were proved to be accurate as well. In a series of works done about ten 



years ago we claimed an accuracy of 3% for Tl [Q, 8% for Pb and 11% for Bi |jT5|. All these 



PNC effects were recently measured to an accuracy of about 1% and found to be in good 
agreement with our predictions. Tliis means tliat our estimates for the theoretical accuracy 
were correct and probably even too pessimistic. For example, in our first calculation of the 
Fr energy levels |jTO| we claimed the accuracy of our predictions to be about 0.5% while the 
actual agreement with latter measurements was found to be 0.1%. The situation was similar 
for the electromagnetic transitions 6s-6pi/2 and 6s-6p3/2 in Cs (see below). These numerous 
tests give us firm ground to believe that the theoretical error in Ep^c (0) indeed does not 
exceed 1%. 

As can be seen from ([l|) an accurate value of the vector transition polarizability f3 is also 
required for the interpretation of the PNC measurements. There are no direct experimen- 
tal measurements of P and so the value P = 27.0(2)ag calculated in [jl6| was used for the 



interpretation of the PNC measurements. The theoretical ratio of the scalar transition polar- 
izability a = — 268(3)aQ to P {(y/P)theory = —9.93(14) [|1^ was in good agreement with the 
corresponding experimental value {a/P)exper = —9.9(1) [T^ available at that time. Since then 
the ratio (a//3) was remeasured to a very high accuracy: (a//3) = —9.905(11) [0. There 
have also been new, very precise measurements of the lifetimes of the 6pi/2 and 6^3/2 states 
of Cs |T^. This allows us to improve the accuracy in the determination of /3, and thus in the 
interpretation of the PNC measurements, by incorporating the new experimental results into 
our calculations. 

The calculations were done using direct summation over the exact intermediate states 

1 1 

P = — ^[(7g|r|npi/2)(npi/2|r|6g)( _ _ p ) (3) 

-(7g|r|np3/2)(np3/2|r|6s)( _^ _^ )] 

Here {s\r\np) is an effective radial integral for electromagnetic transitions between exact 
atomic eigenstates, which are related to the reduced matrix elements by 



{s\\r\\Pi/2) = {Pi/2\\r\\s) = ^{s\r\pi/2) 



(4) 



/4 

{s\\r\\p3/2) = -{P3/2\\r\\s) = y-(s|r|p3/2) (5) 

It easy to see that P vanishes in the absence of the spin-orbit interaction, which splits energy 
levels and radial integrals. Thus, it is practically impossible to do accurate calculations of 
P using experimental results due to the strong cancelation between different terms, which 
causes the relative statistical error to be larger. Therefore, we calculated the scalar transition 
polarizability a instead and used the measured ratio a/P to find p. Note however that the 
calculation of a/P using theoretical radial integrals and experimental energies reproduces the 
experimental value for this ratio to an accuracy of about 1%. 
The expression for a is given by 

11 

" = ^[(73|r|r2pi/2)(npi/2|r|6g)( _ + _ ) (6) 

+2{7s\r\np3/2){np3/2\r\6s){ _^ + _^ )] 

^7s ^np^/2 ^np^/2 



Here all of the major terms produce positive contributions. This reduces the error in the final 
result. 98% of the value of a is given by the intermediate 6p and 7p states. The 6p state 



practically does not contribute to the error in the final result. Our calculations of the 6s - 
6p electromagnetic amplitudes were recently confirmed to an accuracy of about 0.1% by very 
accurate experimental measurements |jT9|. The 6p - 7s amplitudes are also known from 
with an accuracy of 0.5% and they agree with the theory. 

The main source of error is the contribution of the 7p intermediate state. The radial inte- 
grals (6s|r|7pi/2) and (6s |r| 7^3/2) are anomalously small due to cancelations between different 
areas of the integration in the single-particle amplitudes. These cancelations substantially 
increase the relative error in the calculated results. Because of this we use the experimental 
values of the 6s - 7p transition amplitudes, which have an accuracy of about 0.7% In 
the relative oscillator strengths were measured using the lifetime of the 6pi/2 state measured 
in [22 1 as a normalization point. Recent measurements of the lifetime are more accurate |19 . 
Therefore we rescaled the experimental 6s - 7p amplitudes from |^ using the new normal- 
ization. Note that the difference between (6s|r|7pi/2) and (6s|r|7p3/2) can be calculated very 
accurately. This is because it is proportional to the mixing between 7p and 6p states by the 
spin-orbit interaction. Indeed, perturbation theory in the spin-orbit interaction ^ gives 

(6s|r|7pi/2) - (6s|r|7p3/2) ~ ^7pep (^Qs\r\6p) H (7) 

^7p ~ ^6p 

The values of the energy levels and spin-orbit splitting can be reproduced almost exactly 
in the numerical calculations by introducing factors into the correlation potential S (since 
the accuracy of the ab initio calculations is high, these factors are close to 1 anyway). The 
calculated matrix element (6s|r|6p) practically coincides with the value obtained from the 
accurate measurements of Ref. Therefore, we believe that the absolute accuracy in the 
calculation of the difference between the doublet radial integrals is always higher than the 
experimental accuracy (to avoid confusion we should note that we use Dirac wave functions, 
i.e. we do not expand in ^ while doing calculations). Thus we can take the experimental value 
of (6s|r|7pi/2) which is measured more accurately, and find (6s|r|7p3/2) using the calculated 
difference {6s\r\7p3/2) — (6s|r|7pi/2). Surprisingly, the result of this procedure gives precisely 
the result of the measurement of the (6s|r|7p3/2) amplitude, which formally has a larger error 
(1.8%). The ratio (6s|r|7p3/2)/(6s|r|7pi/2) also has a smaller experimental error (0.4%) than 
the error in (6s|r|7p3/2) ||21|. Therefore, we may assume that the actual relative error in the 
(6s|r|7p3/2) is 0.7%, similar to that in (6s|r|7pi/2). We use theoretical values of the (7s|r|7p) 
transition amplitudes since we believe that the expected theoretical error here ( 0.3%) is 
smaller than the experimental error. All higher transitions, including continuum and core 
electron transitions, were also calculated theoretically, even though their contribution was 
small (see below). 

The result of the calculation of a is as follows 

a = a(6pi/2) + 0(6^3/2) + a{7pi/2) + 0(7^3/2) + a{others) = (8) 
-32.39(0.17) - 92.56(0.46) - (37.79 + 103.01)(1.14) - 3.25(0.20) = -269.0(1.3) 

We used experimental energy levels from and radial integrals from Table | to calculate 
the contributions of the 6p and 7p states. We used both experimental and theoretical data 
to select the "best values" of these integrals. Note that the errors in a(7pi/2) and a(7p3/2) 
are proportional and so we added them. When new data for electromagnetic amplitudes are 
available it will be easy to refine this result by multiplying the corresponding term by the 
ratio of the new amplitude to the old one. 

This value of a combined with the measurements oi a/(3 gives 



p = 27.15(13)a|^. (9) 



The result of the direct calculation using radial integrals from table | is /3 = 27.00. The 
results of other works are (3 = 27.0(2) H, /3 = 27.2(4) [||], (3 = 27.3(4) H, /3 = 27.17(35) 



3^ . Using (j^) , the measurement (0) , the mean value of the theoretical amplitudes and 



|e|/aQ = 5.1422 x lO^'^mV/cm, we obtain 

Qwiexper) = -72.4:l{25)exper{80)theor- (10) 

Comparing this result for with the theoretical value p7 



Qw{theor) = -73.20(13) - 0.85 - 0.005T, (11) 

we can find the Peskin-Takeuchi parameter S characterizing new physics beyond the Standard 
model (i.e. weak isospin conserving radiative corrections produced by new particles) 

S + 0.006T = -1.0(0.3)exper(1.0)tW (12) 



We can also use the calculation of the extra Za;-boson contribution in the 5*0(10) model |27] 



AQw = 0.4(2iV + Z){^r = 8AA{^f (13) 

Mz, Mz, 

to find the limit for the mass of this boson 

Mz^ > 550GeV (14) 

The natural question is: can we refine the value of the Ep^c calculation using experimental 
ii^l-amplitudes? Unfortunately, the experimental accuracy at the moment is not good enough 



to give an improvement. For example, we can use the results of the work [|T6|, where the 
direct sum-over-states approach was discussed in detail. The theoretical result of the direct 
summation was 

EpNc = -0.907(9)10-i^z|e|ao(-^) (15) 



Replacing the ii^l-amplitudes calculated in [T^ (see table IV of that work) with the values 
from Table | gives 

EpNC = -0.902(ll)^i(~ 7)oi;,e40-^^2|e|ao(-^) (16) 

Here we separated the error coming from the 6p and 7p ii^l-amplitudes from the error coming 
from all other sources, including the weak matrix elements and the amplitudes for transitions 
to the states above 7p. The error in the weak matrix elements can be roughly estimated using 
the deviation of the calculated hyperfine intervals from the experimental values since both the 
weak and hyperfine interactions are approximately proportional to the density of the electron 
wave function near the nucleus. Note that, the error from the i?l-amplitudes exceeds the 
error in the theoretical values for the Ep^c (0)- To avoid confusion we should stress that 
the calculation in Ref. was based on the Green's function technique and does not contain 
partial cancelations of the different terms which increase the error in the direct sum-over-state 
approach. 

In conclusion we would like to stress that accurate measurements of the i?l-amplitudes 
(Table |p are very desirable for an improvement of the interpretation of the PNC measurements 
in Cs. For a the most important improvement would be a more accurate value of the 6s - 7p 
amplitude. An improvement for the 7s - 7p amplitude is also very important because of the 
disagreement between theory and existing data. 

The authors are grateful to David DeMille for helpful comments. This work was supported 
by the Australian Research Council. 
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TABLES 



TABLE I. Radial integrals used in the calculation of a 



np 


(6s r np) 


(np r 7s) 






5 ^c|Cl(')7') 


6P3/2 


-5.4824(62) 


5.605(27) 


A(6p3/2 - 6P1/2) 


0.0267 


0.4154 


7pi/2 


-0.3460(26) 


-12.597(38) 




-0.5040(38) 


-12.372(37) 


A (7^3/2 - 7^1/2) 


-0.158 


0.225 



